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We study the ground-state properties of an extended periodic Anderson model to understand the 
role of Hund’s coupling between localized and itinerant electrons using the density-matrix renor¬ 
malization group algorithm. By calculating the von Neumann entropies we show that two phase 
transitions occur and two new phases appear as the hybridization is increased in the symmetric 
half-filled case due to the competition between Kondo-effect and Hund’s coupling. In the inter¬ 
mediate phase, which is bounded by two critical points, we found a dimerized ground state, while 
in the other spatially homogeneous phases the ground state is Haldane-like and Kondo-singlet-like, 
respectively. We also determine the entanglement spectrum and the entanglement diagram of the 
system by calculating the mutual information thereby clarifying the structure of each phase. 
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I. INTRODUCTION 

The colossal magnetoresistance observed in certain 
manganites has attracted much attention in condensed 
matter physics— *2 Typically perovskite manganites, for 
example I^/aCa^MnOa, exhibit this enormous en¬ 
hancement of the magnetoresistance— It is now gener¬ 
ally accepted that the colossal magnetoresistance orig¬ 
inates from the strong correlation between electrons, 
namely, the double exchange and superexchange inter¬ 
actions, Jahn-Teller effect etc. play an important role, 
although this phenomenon is far from being understood.' 5 
In these compounds due to the crystal field, the five-fold 
degenerate d orbitals are split into two-fold degenerate 
( e g ) and three-fold degenerate ( t 2 g ) orbitals which are 
coupled to each other via Hund’s coupling. The elec¬ 
trons in the e g orbitals are delocalized while the t 2 g elec¬ 
trons are localized.— To understand the compounds in 
question, it is crucial to investigate models with more 
than one orbital. Hybridization between different or¬ 
bitals is claimed to be important also in other com¬ 
pounds like Ca 2 _ a; Sr a; Ru 04 — This compound exhibits 
several remarkable phenomena like orbital-selective Mott 
transition^ or heavy-fermion^ behavior as the chemical 
concentration or temperature is varied. To account for 
these effects, multiorbital systems have been in the focus 
of intensive research recently— r— Our purpose here is to 
consider further the properties of multiorbital systems. 

The starting point of our investigations is the periodic 
Anderson modelji^r— which is a minimal model with two 
kinds of electrons: 
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Here ^k a (c fccr ) creates (annihilates) an electron with spin 
a in a wide band with dispersion curve £k■ This band 
will be described in the tight-binding approximation with 
nearest-neighbor overlap only. Furthermore, (fj a ) 
creates (annihilates) localized electrons at site j with 
spin g and the corresponding particle number operator 
is h^J = fj a fj a ■ They can be mixed to the states of the 
wide band by the matrix element 14. When the states 
of the wide band are written in a real-space represen¬ 
tation, this mixing is in general non-local, however, we 
neglect its fc-dependence in the following, 14 = V. The 
Coulomb interaction between localized electrons is Uf, 
and their on-site energy is £/. 

Our goal is to study the interactions between the two 
kinds of electrons, whose general form in momentum 
space can be written as: 

U = HpAM + ^7 X U ^'^ C k+qafl-qa'fk'a' C ka^ ( 2 ) 

kk' q 
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which is visibly non-local in real-space for a general 
Ucrcr'iq) interaction. In what follows we consider only 
on-site interactions between the two kinds of electrons, 
taking into account the interorbital Coulomb interaction 
and the direct exchange between them, described by the 
rotationally invariant Kanamori Hamiltonian—* In addi¬ 
tion to that we also include an on-site repulsion for the 
delocalized electrons. To write the Hamiltonian in a con¬ 
venient form, we use the real-space representation for 
the creation and annihiliation operators of the electrons 
in the wide band (cj CT , c^), then the total Hamiltonian 
becomes: 

H = Hpam + U c J2 + X ( U °f ^ 

3 j™' 
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where U c is the Coulomb repulsion within the wide band, 
U c f denotes the local interaction between the two kinds of 
electrons and J is the Hund’s coupling. It is worth noting 
that one can always diagonalize the bilinear part of the 
Hamiltonian obtaining two orthogonal bands. However, 
in that case the intraorbital interaction becomes non¬ 
local in terms of the original operators— 

The role of the term U c and U c f has been examined 
in several papers— - — It has been shown, that U c leads 
to significant enhancement of the effective mass in the 
Hondo regime, while U c f causes critical valence fluctua¬ 
tions in the mixed valence regime and can lead to charge 
ordering in infinite spatial dimensions. The full Hamil¬ 
tonian in Eq. m including J has also been investigated 
thoroughly by dynamical mean-field theory (DMFT) in 
infinite spatial dimensions,— ' i 16 i 26 i 27 and it was revealed 
how the Hondo and Mott insulating states compete with 
the metallic state in the half-filled case if the system is 
assumed to be paramagnetic— £ It turned out, however, 
that magnetic long-range order can also be present in 
the model, namely two types of antiferromagnetic or¬ 
der emerge beside the Hondo insulating state— - The oc¬ 
curence of these phases originates from the competition 
between Hund’s coupling and Hondo effect. While the 
former one aligns the spins of localized and itinerant elec¬ 
trons ferromagnetically at a given site, the latter one tries 
to screen the localized spins by forming singlets with the 
itinerant electrons. It has also been discussed how the 
hybridization affects the orbital-selective Mott localiza¬ 
tion emerging in two-orbital Hubbard-models— - 

Since the DMFT approach completely neglects the spa¬ 
tial fluctuations, which is only valid in infinite dimen¬ 
sions, it is necessary to investigate low-dimensional sys¬ 
tems where quantum fluctuations are known to be much 
stronger. Our main goal in this paper is to explore the 
one-dimensional behavior of the Hamiltonian in Eq. ©. 
In earlier papers it has been show n 25 i 28 i 29 that there is 
no quantum phase transition in one dimension in the 
absence of the Hund’s coupling in contrast to the in¬ 
finite dimensional case. The competition between the 
Hund’s coupling and the Hondo effect may lead to the ap¬ 
pearance of quantum phase transitions and unexpected 
phases even if true long-range order cannot be present in 
one dimension. 

We apply the density-matrix renormalization-group 
method^ - — (DMRG), which is a powerful tool to find 
the ground state and to determine the correlation func¬ 
tions. Further advantage of the DMRG method is that 
we can easily determine the von Neumann entropies^ - — 
of single and multisite subsystems, without the need to 
calculate excited states, which is in general difficult near 
a critical point, and their anomalies can be used to detect 
quantum phase transitions^ 1- — 

In our DMRG calculation we applied the dynamic 
block-state selection algorith m 45 ' 46 in which the thresh¬ 
old value of the quantum information loss, y, is set a 
priori. We have taken y = 3 • 10 -6 . A maximum of 
2000 block states is needed to achieve this accuracy, and 


the truncation error was in the order of 10~ 7 . Such low 
value of y is necessary in order to obtain ’smooth’ data 
sets close to critical points. We investigated chains up to 
a maximum length L = 120 with open boundary condi¬ 
tions and performed 8 sweeps. 

The setup of the paper is as follows. In Sec. II. we 
define the von Neumann entropies of various subsystems 
used in our analysis, and the mutual information . 38 i 47 ' 48 
In Sec. III. A, B and C we discuss the properties of the 
phases occuring in the model using the mutual informa¬ 
tion and the eigenvalue spectra of the two-site density 
matrices. In Sec. III. D we discuss the differences be¬ 
tween the phase diagram obtained in the DMFT and for 
the one-dimensional model. Finally, in Sec. IV. our con¬ 
clusions are presented. 

II. VON NEUMANN ENTROPIES 

The von Neumann entropies of different subsystems are 
known to exhibit anomalies near critical points i 39 ' 41 ' 42 
We examined the one-site s^, two-site entropies and 
the block entropy which is the entropy of the subsys¬ 
tem containing sites from 1 to L/2. These quantities 
can be obtained from the appropriate reduced density 
matrices—i 36-39 The entropy of a single site can be ob¬ 
tained as 

Si = -Trpilnpj, (4) 

where pi is the reduced density matrix of site i, which 
is derived from the density matrix of the total system 
by tracing out the configurations of all other sites. We 
also define the entropies corresponding to the two types 
of electrons at a site (s\ c \ s^) in the following way: 

s\ c) =-Trp^lnp^, (5) 

s\ f) =-Trp ( j f hnp ( /\ ( 6 ) 

where pf' 1 {p\^) is obtained by performing an additional 
trace over the remaining / (c) degrees of freedom at site 
i. The two-site entropy is written as 

Sij — Tr Pij In Pi] • (f ) 

where p^ is the two-site reduced density matrix of sites 
i and j. We can also introduce the partial two-site en¬ 
tropies for type a electrons on site i and type b electrons 
on site j: 

s (ah) = _ Trp (a 6 ) ln p (ab), ^ 5 £ { c , /} ( 8 ) 

where p is derived from p^j by tracing out the states of 
the other electrons. The mutual information^ - — which 
measures the entanglement between sites i and j can be 
obtained from: 

Iij — Sj V Sj Sij , 


( 9 ) 
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while the mutual information between a and b type elec¬ 
trons on sites i and j is defined as 


A ab ) — «,(“) 

b i 


Jb) 


„(° b ) 


( 10 ) 


which measures all correlations both of classical and 
quantum origin between a and b type electrons on sites 
i and j. In what follows we refer to as the entan¬ 
glement between these components. Finally, the block 
entropy is defined as 


s(L/2) = -Trp A \np A , (11) 

where A denotes the subsystem which contains the sites 
from 1 to L/2. In contrast to the one- or two-site 
entropies, which have a finite upper bound, the block 
entropy grows as O(lnL) for one-dimensional critical 
systems] 36 ’ 37 



FIG. 1. The block entropy, s(L/ 2) as a function of hybridiza¬ 
tion for different chain lengths and J/U = 0.1. The lines are 
guides to the eye. 


III. RESULTS 

In what follows we consider the half-filled case and 
nearest-neighbor hopping between delocalized electrons, 
e(k) = —2fcosfc, and use the half bandwidth, W = 21, 
as the energy scale of the system and set £f = 0. For 
simplicity we assume U c = U / = U and U = U c f + 2 J. 
In the absence of Hund’s coupling the ground state in one 
dimension is either a collective singlet or consists of less 
entangled local Kondo-singlets depending on the values 
of Coulomb interactions and the hybridization^ There is 
no quantum phase transition between these phases, just a 
smooth crossover separates them. To examine the effect 
of the Hund’s coupling, firstly we consider what happens 
for a finite Hund’s coupling, namely for J/U = 0.1 and 
0.3, with U = 4 W as the hybridization is varied. 

Firstly, we investigate the block entropy of one half of 
the chain. This quantity is a smooth function of V for 
any U when J = 0. For any finite J, however, two peaks 
appear in the block entropy as can be seen in Fig. |T] for 
different chain lengths for a fixed value of J/U = 0.1, 
where the two peaks are around V/W = 0.57. It is 
clearly observed that the height of the peaks increases as 
the chain length is increased. We know that maxima in 
the block entropy can be attributed to quantum critical 
points^ if they evolve into anomalies in the thermody¬ 
namic limit. Two peaks may indicate the existence of 
two phase transitions separating three different phases. 
To check if it is indeed the case, one has to show that the 
peaks remain separated and do not merge in the thermo¬ 
dynamic limit. The finite-size scaling of the position of 
the peaks is shown in Fig. [2] We could treat systems 
with L = 60 sites near the critical points due to the high 
value of the block entropy. To determine the positions of 
the maxima accurately we used a cubic spline interpola¬ 
tion. Figure [2] shows that the position of the peaks as a 
function of 1 /L can be fitted well with a linear function 
and the phases remain separated in the thermodynamic 



FIG. 2. The main figure shows the finite-size scaling of the 
peaks occurring in the block entropy for J/U = 0.1, while the 
inset shows the same for J/U = 0.3. The lines are the best 
linear fits to the data. The estimated error is comparable to 
the size of the symbols. 


limit having positions: 

Vr/W = 0.5677(0) 

V^/W = 0.5704(1) 1 ’ 

for J/U = 0.1. As a remark we mention that to estimate 
the error of the data points in Fig. [3] we use the follow¬ 
ing analysis. The quantum information loss was a priori 
set to x = 3 • 10~ 6 . Therefore the relative error of the 
ground state energy is estimated to be of the same order 
of magnitude, 5E re \ ~ 10 -6 . Since the ground state ener¬ 
gies are of the order of E ~ 10 2 , the absolute error of the 
energies is of the order of <5-E a bs ~ 10 -4 , so the error of 
the block entropy can be expected to be not larger than 
bSabs ~ 10 -4 . The block entropy values near the critical 
points were sampled with an equidistant step 5-10~ 4 (for 
better visibility not all are shown in Fig. [T]). The error of 
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the spline fit is also expected to be 10 -4 . All in all the 
overall error is expected to be 10 ~ 4 , which is comparable 
to the size of the symbols in Fig. [2] 

We repeated the same calculation for other values of 
J/U. Our results indicate that two separate peaks are 
present in the block entropy for any finite J. The distance 
of the peaks is shown in Table Q] for two values of J/U 
for different chain lengths. It is clear that their distance 


J/U 

L = 16 

L = 24 

L —y 00 

0.1 

0.016 

0.011 

0.0027(1) 

0.3 

0.040 

0.028 

0.006(8) 


TABLE I. The distance of the peaks of the block entropy for 
several chain lengths and J/U ratios. The extrapolation was 
performed using chains up to L = 60. 

grows as J is increased, furthermore, the critical values of 
V are also shifted towards larger values. In what follows 
we analyze the ground state properties of each phase for 
J/U = 0.1 using the mutual information and correlation 
functions to check if the peaks indeed separate different 
phases. 

A. The Kondo singlet phase for Vf 1 < V 

Firstly, we consider what happens for large hybridiza¬ 
tion where the effect of Hund’s coupling is expected to be 
small compared to that of hybridization and the Coulomb 
interaction and the properties of the J = 0 model are ex¬ 
pected to be recovered. We examine how the individual 
system components are entangled to each other using the 
mutual information. We have seen already in Fig. [T| that 
for large values of the hybridization the block entropy 
decreases rapidly, which indicates a less entangled state. 
This is the case indeed, as is seen in the entanglement 
map in Fig. [3] We can see that very strong on-site entan¬ 
glement appear, while the one-particle states on different 
sites are hardly entangled to each other. To describe 
the physical origin of this difference in the entanglement 
within a site and between neighboring sites we calculated 
the eigenvalues (w 7 , 7 = 1 ,..., 16) of the two-site den¬ 
sity matrix pf^ and the corresponding eigenfunctions. 

From p ^/2 l /2 we f° un( i that for V/W = 0.8 one of its 
eigenvalues is larger by two orders of magnitude than the 
others and the corresponding eigenvector is: 

^l/Il/2 = 0-5574(| t)g 2 | | +>g 2 | t )%) 

+ 0.4350(| n>g a |0)i ? 2 + |0)g a | U)l/ 2 )- 

Here |0)“, | f)“> I 4-)“) I t4-)“ denote the four possible 
states of electron type a on site i. We can see that strong 
on-site singlets are formed between localized and delocal¬ 
ized electrons, which we may refer to as Kondo-singlets, 
since this is the consequence of the enhanced Kondo ef¬ 
fect. The entanglement between nearest-neighbor sites is 



FIG. 3. Schematic view of all components of the mutual in¬ 
formation (4f\ 4 c/) , 4/ J) ) for V/W = 0.8 and L = 16. 
The numbers are the site indices. The inner and outer circles 
denote type / (localized) and type c (itinerant) electrons. 

smaller by two orders of magnitude than the on-site en¬ 
tanglement between localized and delocalized electrons. 
Therefore, the ground state is almost a product state. 
Since the eigensystem of l/2+i an< ^ P'l/2. L/2+ 1 4S 
quantitatively the same, we consider only the former 
one. The eigenfunction belonging to the most significant 
eigenvalue of pf^l £/ 2 +i rea ds: 

Tec) 

< Pl/ 2 ,l/ 2 +i ~ 

0.6900(| t)g 2 | l)g 2+1 - I l)g 2 | t}g 2+1 ) (14) 

+ 0.1545(| U)g 2 |0)g 2+1 + |0)g 2 | tl>i% +1 ), 

that is, the nearest neighbor coupling between the spins 
is antiferromagnetic. We checked that the mutual infor¬ 
mation components have their bulk values at L = 16 al¬ 
ready, which is due to the hardly entangled ground state. 
Indeed, the properties of this phase agree with the known 
behavior of the conventional periodic Anderson model for 
large hybridization. 

B. The Haldane-like phase for V < VT 

A new phase is expected to appear for small hybridiza¬ 
tion, where J dominates. Here we discuss the properties 
of the phase emerging for V < V/ T using the mutual 
information. The entanglement diagram containing all 
types of the mutual information is shown in Fig. [I] for 
V/W = 0.3. One can see that the strongest entanglement 
is developed between neighboring delocalized electrons 
and moderately strong entanglement is present between 
more distant sites. 

Firstly, we consider how the eigenvalue spectrum looks 
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FIG. 4. The same as in Fig. [3] but for V/W = 0.3. 

like for L ^ 2 and L = 16. We found that one of its 
eigenvalues is threefold degenerate and larger by an order 
of magnitude than the others. Its value is very close to 
1/3 and the three corresponding eigenfunctions, (jr^j 2 '’2/2 
7 = 1,2,3, read: 

j/c /). 1 _ I | t \(/) 

™L/2,L/2 — I UL/2\ I/i/2’ 

^L/2,L/2 = y|(l t)i/ 2 l Dl/2 + I 'I')L/21 t)i^/ 2 ))> (1^) 

^L/2,’l/2 ~ I 'I') L/2 1 J-)l/2- 

That is, the electrons on the same site are in a state where 
the S = 1 triplet components have the largest weights. 

As a next step we examine the entanglement be¬ 
tween nearest neighbor sites. Since the eigensystems 

of Pl/1,l/2+v Pl/Il/ 2 +i and Pl/Il/ 2 +i are quantita- 

( f f) 

tively very similar, we only present results for P ^/2 l/ 2 +v 
The eigenfunction corresponding to the most significant 
eigenvalue is: 

A ff ) 

"L/2,L/2+l 

0.7071(1 1)^/2 1 i) l /} 2+1 - I 4 -) L/2 1 t)i7a +1 ) (16) 

+ 0.0014(| U)l/ 2 \°)l/2+i + 10)^1 tDi^+J, 

which means that the entanglement between the neigh¬ 
boring sites results mainly from the singlet formation. 
Furthermore, we investigated how the mutual informa¬ 
tion components scale as the system size is increased. 
This is shown in Fig. [5] It is obviously seen that the 
bonds hardly change as the chain length becomes larger. 
The extrapolation was performed using a quadratic poly¬ 
nomial: 

I% b) (L)=I% b) +A/L + B/L 2 , (17) 

where I^j b \ A and B are free parameters. 


0.4 
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FIG. 5. The finite-size scaling of the mutual information for 
V = 0.3W. The dotted lines are quadratic polynomial fits to 
the data. The inset drawing denotes a segment of the middle 
of the chain, the circles denote the type c and type / electrons 
on a given site, while the symbols denote the bonds shown in 
the main plot. 



This confirms that the spins are aligned ferromagneti- 
cally within a site, but they couple antiferromagnetically 
between nearest-neighbor sites. The former one is a con¬ 
sequence of the strong Hund’s coupling which prefers par¬ 
allel alignment of the spins, while the latter one is due to 
the RKKY-interaction mediated by the conduction elec¬ 
trons. The same structure was confirmed for V = 0. 
These findings suggest that the model in this regime can 
be considered as an S' = 1 Heisenberg chain with an¬ 
tiferromagnetic nearest-neighbor coupling, therefore this 
is a Haldane-like phase. The ground state is a singlet, 
however, in the thermodynamic limit the ground state of 
an open chain becomes degenerate with the first S = 1 
excited state due to the end spins. This is a well-known 
property of the Haldane phased 


C. The dimerized phase for Vi T < V < V™ 

Finally, we examine the properties of the narrow in¬ 
termediate phase, whose appearance is indicated by the 
analysis of the block entropy. Using the tools applied in 
the previous subsections we examine the spatial struc¬ 
ture of the ground state. The entanglement map shown 
in Fig. 1 is drastically different from that in the previ¬ 
ous phase. It is clearly seen that strong and weak bonds 
alternate along the chain, which suggests a spatially inho¬ 
mogeneous, dimerized ground state. Before investigating 
the physical processes that contribute to the creation of 
the strong entanglement, it is important to check if the 
dimerization remains finite in the thermodynamic limit. 
We can introduce two types of order parameters for the 
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FIG. 6. The same as in Fig. |3] but for V/W = 0.569. 

dimerization: 

D[ a) (V) = lim 

L—tcx 

D 2 (V)= lim |s(L/2, V) — s(L/2 + 1, V)\ . (19) 

L—> oo 

Since D[ a \V) is a local quantity, we expect that it is 
less sensitive to the boundary effects. It requires, how¬ 
ever, the calculation of several correlation functions, and 
their computation time scales as L 2 which can be com¬ 
putationally demanding. The quantity in Eq. m is 
computationally less demanding, but since the block en¬ 
tropy is a non-local quantity its convergence to the bulk 
value may be slower. Instead of showing (E) directly, 

we plot the individual values of the mutual information 
components ( 4 / 2 ,l /2+1 , 4 “-i,i/ 2 ) and investigate their 
size-dependence. In this case we could consider chains up 
to L = 120, since in the intermediate phase the block en¬ 
tropy has a much lower value than near the critical points 
and its low value also indicates a less entangled ground 
state as expected for a dimerized phase. This is shown 
in Fig. [H In Fig. [7] we used a quadratic polynomial, Eq. 
m for the extrapolation, which gives an upper bound 
for the bond strengths. The data can also be fitted using 
a power-law function: 

l\f\L) = l\f ) +A/L b , (20) 

where , A and B are free parameters. The residual 
sum of squares is roughly of the same order of magni¬ 
tude for both types of fits, ©(lO” 6 ), therefore we give 
the values of the order parameters for both fits in Ta¬ 
ble EH The quadratic extrapolation clearly overestimates 
the order parameters while the polynomial fit underesti¬ 
mates it, since we expect that the order parameter begins 
to saturate as soon as the bulk limit is achieved. 

The other order parameter defined in Eq. m is shown 
in Fig. [8] for J/U = 0.1 and J/U = 0.3. The extrapola- 



typeoffit D[ c) {V) D[ f) (V) P 2 (V) 
quadratic 0.043(0) 0.215(9) 0.32(4) 
power-law 0.03(1) 0.17(8) 0.26(0) 


TABLE II. The extrapolated dimerization order parameters 
for different types of fits at V/W = 0.569. 



l/L 


FIG. 8. Finite-size scaling of the order parameter in Eq. m 
The symbols •, ■ belong to V/W = 0.569, J/U = 0.1 and 
V/W = 1.665, J/U = 0.3, respectively. The dotted and 
dashed lines denote the quadratic and power-law fits, respec¬ 
tively. 


tions using the two different fits are shown in Table EH 
From these calculations we conclude that the intermedi¬ 
ate phase remains dimerized in the thermodynamic limit, 
which is a clear sign of spontaneous symmetry breaking. 
Moreover, it is a surprising phenomenon, since dimeriza¬ 
tion has not been observed at all in the periodic Anderson 
model without the Hund’s coupling. 

Now, we return to the entanglement patterns in Fig. 
E We consider first the entanglement between local- 

( f f) ( f f) 

ized electrons by examining P// 2-1 l /2 an< ^ Pl /2 L/ 2+1 


/(“) (V) — (V) 

1 L/2,L/2+l\ V ) 1 L/2-l,L/2\ V > 


( 18 ) 
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for L = 16. We have seen that due to the finite-size 
effects the strengths of the bonds change, but the qual¬ 
itative picture, which can be obtained from the analysis 
of the density matrices, remains. For pY/ 2-1 l /2 one °f 
the eigenvalues is larger by an order of magnitude than 
the others, and the corresponding eigenfunction is: 

Af f) 

™L/2-l,L/2 — 

0.7071(| tfSJLrl +>% - I +>V2-il t)^) (21) 

+o.oo6(i tDiVii 0 )^ + \°)l/ 2 ~i\ u)%} 2 ), 

which means that the origin of the strong entanglement 
between localized electrons on neighboring sites is the 
singlet formation. If we consider the neighboring bond, 
we obtain from L / 2+1 that one °f the eigenvalues is 

w 1 = 0.3953 and there is a threefold degenerate eigen¬ 
value W 2 = 0.1500. The eigenvector corresponding to the 
former one is essentially the same as in Eq. m, while 
the eigenvectors corresponding to the latter one are the 
triplet components described in (1151) . Due to the fact 
that triplet components are mixed with a larger weight 
to the singlet component, it destroys the singlet bond be¬ 
tween the localized electrons resulting in a much weaker 
entanglement. Qualitatively the above considerations re¬ 
main valid for the explanation of entanglement between 
the itinerant electrons. Lastly we examine the entangle¬ 
ment within a site with the help of p^f 2 L / 2 - I* 1 this case 
we have again a non-degenerate eigenvalue, = 0.2467, 
and a threefold degenerate one, W 2 = 0.2274. The eigen¬ 
vector corresponding to uq is 

2 = 0.5885(| f)i%| - U)i%l t)^) 

+ 0.3921(| U) { l/2\°)l/2 + l°)i C /2l U)%), 

while the eigenvectors of w 2 are the triplet states in (fl5l) . 
It can be seen easily, that the on-site spin correlation is 
still ferromagnetic, but significantly reduced compared to 
the Haldane-like phase. While the on-site singlet state 
has negligible weight in the Haldane-like phase, in the 
dimerized state the on-site triplet and singlet states are 
mixed with comparable weights. 

D. Discussion 

In the light of the above results it is worth examining 
the nature of the phase transitions and comparing the 
properties of the phases to what has been obtained in 
infinite dimensions. 

As we mentioned, there is no phase transition when 
J = 0, where the ground state is Kondo-singlet-like dis¬ 
cussed in Sec. III. A. For any finite J two new phases 
appear, namely, a Haldane-like and a dimerized phase 
whose properties are discussed below, and they disap¬ 
pear as J —> 0. 


We have seen that for V < Vj cr the ground state is 
Haldane-like while for V > V 2 r Kondo-singlet-like, and 
both are gapful and homogeneous. For V 1 cr < V < V/ 1 
the translational symmetry is broken due to the dimer¬ 
ization. It is worth noting that similar phase diagram has 
been obtained in frustrated spin ladders^ where on-site 
and nearest-neighbor antiferromagnetic couplings com¬ 
pete with each other. 

The existence of different phases can be corroborated 
by investigating the entanglement spectrum^ which is 
known to be a useful tool to detect symmetry-protected 
topological order^ It is obtained from the eigenvalues 
(Aj) of the density matrix of a half chain. In the Haldane- 
like phase the degeneracy of each eigenvalue is an even 
number while both even and odd degeneracies may 
appear in other phases. For short chains both even and 
odd degeneracies occur for any V/W. This picture drasti¬ 
cally changes for longer chains, namely, for L > 90. Here 
the correlation between the end spins suddenly vanishes, 
which is marked by a jump in the block entropy for small 
V/W. While even and odd degeneracies remain in the 
spectrum for V > Vj cr , the entanglement spectrum con¬ 
tains nearly degenerate eigenvalues whose multiplicity is 
even for V < V 1 cr . In the thermodynamic limit the eigen¬ 
values becomes exactly degenerate. In Fig. 0 the low- 
lying eigenvalues and their degeneracies are shown with 
a logarithmic scale corresponding to the three different 
phases for L = 120. While the spectra of the dimerized 

Haldane-like Dimerized Kondo-like 




V/W = 0.3 V/W = 0.569 V/W = 0.8 

FIG. 9. The entanglement spectrum in the three different 
phases for J/U = 0.1 and L — 120 (for better visibility only 
the lower part is shown). 

and Kondo-like phases contain eigenvalues with odd and 
even degeneracy, only even degeneracy is present in the 
Haldane-like phase in agreement with our expectation. 

It is interesting to compare these findings with what 
has been obtained by DMFT.— Surprisingly, three dis¬ 
tinct phases appear also in the DMFT phase diagram, 
although their properties are significantly different. In 
DMFT there is a phase with antiferromagnetic long- 
range order in which the on-site spins are coupled ferro- 
magnetically while the nearest-neighbor coupling is anti- 












ferromagnetic (AF I phase). Further increase of the hy¬ 
bridization drives the system into an intermediate phase 
(V 1 cr < V < V£ T ) where another type of antiferromag¬ 
netic order takes place. Here the on-site coupling be¬ 
comes antiferromagnetic, while the nearest-neighbor cou¬ 
pling remains antiferromagnetic (AF II phase). Finally, 
for V£ r < V Kondo-like behavior is realized. 

In one dimension we cannot expect true long-range 
magnetic order, only slow decay of the correlation func¬ 
tions. For V < Vj cr the on-site spins are parallel due 
to the strong Hund’s coupling. This Haldane-like phase 
might be the residue of the AF I ordered phase obtained 
in the DMFT calculation. Above the second critical 
point, V™ < V, the ground state is homogeneous, and 
strong on-site correlations appear, which originate from 
the enhanced Kondo-effect and the sites are occupied 
more and more by two localized or delocalized electrons 
or vice versa. The properties of the Kondo phase are 
consistent with what has been obtained in DMFT. Both 
methods exhibit an intermediate phase between them, 
however, their properties are completely different, which 
is caused by the enhanced quantum fluctuations. The ap¬ 
pearance of the dimerized phase may originate from the 
competition between the Haldane-like and Kondo singlet 
phase. This may not surprise us if we recall that dimer¬ 
ization has been found in a two-component system via a 
purely electronic mechanism^ The present case is sim¬ 
ilar although for antiferromagnetic interactions between 
localized electrons. 


IV. CONCLUSIONS 

We investigated an extended periodic Anderson model 
in one dimension to understand the role of the Hund’s 
coupling between itinerant and localized electrons. We 


carried out accurate DMRG calculations with quantum 
information loss y = 3 • 10 -6 up to the accessible chain 
length L = 120. For such error margin we found that 
the two competing processes, the Hund’s coupling and 
the Kondo effect, lead to the appearance of two new 
nontrivial phases. The enhanced quantum fluctuations 
in one dimension crucially affects the properties of the 
phases obtained in infinite spatial dimensions. We per¬ 
formed a quantum information analysis and determined 
the entanglement spectrum in the various phases and the 
entanglement patterns between the system components. 
Moreover, using the eigensystem of the two-site density 
matrices we examined what physical processes lead to 
the development of entanglement. We found that for 
V < Vf 7 the itinerant and localized electrons form a local 
triplet within a site and couple to the nearest-neighbor 
sites antiferromagnetically. Here the model can be con¬ 
sidered as an S' = 1 spin chain with antiferromagnetic 
coupling and we have a Haldane-like ground state. This 
is also corroborated by the entanglement spectrum. For 
V 1 cr < V < U 2 cr we have an intermediate dimerized phase, 
in which strong and weak singlet bonds alternate between 
nearest neighbor itinerant and localized electrons. We 
have seen that the region where dimerization occurs ex¬ 
pands as the Hund’s coupling is increased and shifts to 
large values of the hybridization. Finally, for V > 
the ground state is homogeneous, and local singlets are 
formed at each site either by two localized/conduction 
electrons or by a localized and an itinerant electron. 
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